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We have carried out a determination of the magnetic-field-temperature (H-T) phase diagram for 
realistic models of the high field superconducting state of tetragonal S^RuCU with fields oriented 
in the basal plane. This is done by a variational solution of the Eilenberger equations. This has 
been carried for spin-triplet gap functions with a d- vector along the c-axis (the chiral p-wave state) 
D ' and with a d-vector that can rotate easily in the basal plane. We find that, using gap functions 

00 , that arise from a combination of nearest and next nearest neighbor interactions, the upper critical 

field can be approximately isotropic as the field is rotated in the basal plane. For the chiral d- 
vector, we find that this theory generically predicts an additional phase transition in the vortex 
state. For a narrow range of parameters, the chiral d-vector gives rise to a tetracritical point in 
, the H-T phase diagram. When this tetracritical point exists, the resulting phase diagram closely 

resembles the experimentally measured phase diagram for which two transitions are only observed 
O ' in the high field regime. For the freely rotating in-plane d-vector, we also find that additional phase 

»' j , transition exists in the vortex phase. However, this phase transition disappears as the in-plane 

f"S !' d-vector becomes weakly pinned along certain directions in the basal plane. 
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I. INTRODUCTION 



It is widely believed that tetragonal Si^RuC^ 1-3 is a spin-triplet chiral p-wave superconductor. In particular, 
a pairing state characterized by a gap function d = z(k x ± ik y ) best explains existing experimental results. The 
observation of the appearance of local magnetic moments below the superconducting transition temperature by muon 
spin-resonance measurements of Luke et al. (fiSR) 4 can qualitatively be accounted for by the two-fold degeneracy of 
the order parameter (a non-degenerate order parameter cannot give rise to local magnetic moments). The nuclear 
magnetic resonance 5 ' 6 and spin polarized neutron scattering measurements 7 carried out for fields applied perpendicular 
s ! i to the four-fold symmetric c-axis show that the spin susceptibility is unchanged by the normal to superconductor 
transition. This is naturally explained by d = z(k x ± ik y ) since d-vector is perpendicular to the magnetic field for 
which no change of spin susceptibility is expected. Also, the Josephson experiments of Liu et al. arguably place 
the strongest constraint on the orientation of d- vector 8,9 and also implies a chiral p-wave superconducting state with 
, d || z. Finally, the observed field distribution of the vortex lattice for the field along the c-axis is not consistent with 
■ a non-degenerate (single component) order parameter but can be accounted for by the chiral p-wave state 10,11 . 

These experiments provide a convincing picture in favor of a chiral p-wave superconductor. However, there are 
some experiments that do not directly support this state. In particular, the more recent Knight shift measurements 
of Murakawa et al. have been carried out for the field along the c-axis 6 . These measurements reveal no change in 
the spin susceptibility. For this field orientation, this would lead to the conclusion that d-vector is in the basal plane, 
not along the c-axis as would be the case for the chiral p-wave state. The simplest interpretation of this experiment 
^ , is that the magnetic field is sufficiently strong as to rotate d from z to the basal plane. This would imply that the 
transition temperatures for d in the plane are close but slightly less than that for d along z. This is possible if 
spin-orbit coupling is weak. Another explanation for the Knight shift data is that the d-vector is in the basal plane, 
O ■ but free to rotate in the plane. This would require all four possible in-plane degrees of freedom to be degenerate (or 
at least nearly degenerate) . If this is the case then for any in-plane field orientation, d- vector will have an in-plane 
• , component perpendicular to the field. Consequently, the spin susceptibility will remain unchanged for fields applied 
y\ 1 in the basal plane as well. 

Another difficulty with the chiral p-wave state is that, for magnetic fields applied in the basal plane, there are two 
qualitative predictions for which there is little experimental evidence. These are: 

(i) the existence of an anisotropy in the upper critical field as the field is rotated perpendicular to the four-fold 
symmetry axis that does not vanish asT-> T c 12 ; 

(ii) the existence of a phase transition in the vortex state in addition to the usual transitions at H C 2 and H c \. This 
additional transition is due to a change in the structure of the order parameter 13 . 

The primary goal of this work is to understand if there are microscopic theories of the chiral p-wave state that 
can lead to situations where the above (i), (ii) predictions do not hold. We find that it is plausible that one of the 
two predictions does not hold, but less likely that both do not hold. Intriguingly, this analysis also points to the 
possibility of a tetra-critical point in the H-T phase diagram. This tetracritical point has features that agree with 
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recent experimental measurements in high magnetic fields 14 . Given the difficulties that the chiral p-wave state has 
explaining the H-T phase diagram, we also address the possibility of an in-plane d-vector to see if it can account 
for the observed phase diagram. This follows a discussion of the limited conditions for which an in-plane d-vector is 
consistent with experimental results. We find that an in-plane d-vector that is nearly free to rotate in the plane can 
explain the H-T phase diagram. 

The paper begins with an overview of Ginzburg Landau theory for the chiral p-wave state to provide the origin of 
the two predictions (i) and (ii) above. Then we discus the role of spin-orbit coupling on the orientation of the d-vector. 
This discussion motivates an examination of the the quasi-classical equations in a magnetic field for which more than 
one irreducible representation is important. The quasi-classical equations are solved using an approximation that is 
valid in high field regime for all temperatures. We present the resulting H-T phase diagrams for the chiral p-wave 
state. Finally, after a discussion of the consistency of an in-plane d-vector with existing experimental results, we 
present results on the H-T phase diagram for this case as well. 

II. GINZBURG LANDAU THEORY 

The simplest framework within which the role of magnetic fields on the chiral p-wave state can be understood is 
the Ginzburg Landau theory. Here we give a brief overview of this theory to demonstrate the origin of the additional 
transition and the anisotropy in the upper critical field that we will discuss later within a microscopic theory. The free 
energy density for the E u representation of D 4 h with a basis r\ = (t] x ,r] y ) [this basis has the same rotation properties 
as (x,y)] is given by 12 ' 15 

/ = -H 2 + |r?| 4 /2 + foiVxVl - VyV* x ) 2 /2 + fo\v*\ 2 \Vy\ 2 + |A^x| 2 + \D yVy \ 2 

+K 2 (\D y r lx \ 2 + \D x r) v \ 2 ) + k 5 (\D zVx \ 2 + |L> Z %| 2 ) (1) 
+K 3 [{D x r] x )(D y T] y )* + h.c] + K 4 [(D y r) x ){D x r] y )* + h.c] + h 2 /(8n), 

where Dj = Vj — j^-Aj, h = V x A, and A is the vector potential. There are three possible homogeneous phases 12,15 : 

(a) 77 - (M)/\/2 (ft > and ft > ft/2), (b) V = (1,0) (ft > and ft < ft/2), and (c) V = (1, l)/y/2 (ft < and 
ft < 0). Phase (a) is nodeless (if the Fermi surface has the same topology as a cylinder) and phases (b) and (c) have 
line nodes. Weak coupling theories give rise to phase (a): the chiral p-wave phase. The application of a magnetic field 
in the basal plane breaks the degeneracy of the two components r\ x and r\ v . For the chiral p-wave state, symmetry 
arguments imply that the vortex lattice phase diagram contains at least two vortex lattice phases for magnetic fields 
applied along any of the four two-fold symmetry axes: {(1, 0, 0), (0, 1, 0), (1, 1, 0), (1, — 1, 0)} 13 ' 16 . To illustrate the 
origin of these phase transitions, consider a zero-field ground state rj = and a magnetic field applied along 

the (1,0,0) direction. Due to the broken tetragonal symmetry, the degeneracy of the rj = (1,0) and the rj = (0, 1) 
solutions is removed by the magnetic field. Consequently, only one of these two possibilities will order at the upper 
critical field. However, if the system is spatially uniform along the magnetic field, then the solution near the upper 
critical field will exhibit a symmetry that the zero-field solution does not. For our example, this symmetry is cither 
a x (if f] = (0, 1) orders at H c2 ) or — a x = U(ir)<j x (if r] = (1,0) orders at H c2 ) where U(ir) is a gauge transformation 
and o x is a reflection through the x-axis. The only way this can occur is if there is an additional phase transition 
as magnetic field is reduced to break this symmetry. The only difference that occurs for the field applied along the 
(1,1,0) direction is that the solution near the upper critical field will be either r] = (1,1) or 77 = (1,-1). 

Another result of the E u theory that follows from Eq. 1, originally shown by Gor'kov, is that the upper critical 
field is anisotropic near T c 12 . Such an anisotropy, for which dH c2 /dT\T=T c is not equal for (1,0,0) and the (1,1,0) 
directions, cannot occur for order parameters that have only one complex degree of freedom 12 . The anisotropy 
in upper critical field near T c has been calculated from microscopic calculations for the in-plane fields, along the 
(1,0,0) and (1,1,0) directions, for a gap function of the form d(k) = z[rj x f x (k) + r] y f y (k)] 17 . These calculations 
show that anisotropy is generally much larger than that experimentally observed. However, under certain special 
circumstances, this anisotropy can be small 17 . To examine the lack of anisotropy for the whole temperature range 
requires a microscopic model that goes beyond the Ginzburg Landau theory as is done below. Note that previous 
microscopic studies of the chiral p- wave state for in-plane magnetic fields 18 ' 19 did not show reveal the physics discussed 
here. This was because the order parameter in these works was fixed to have the form 77 = (1, i) for all magnetic fields 
and temperatures. Such an approximation is valid only for fields much smaller than H c2 . 
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III. SPIN-ORBIT COUPLING AND THE ORIENTATION OF d 



An important interaction in determining the specific spin-triplet pairing state in Sr 2 Ru04 is spin-orbit coupling. 
We quantify this in this section. From a symmetry point of view, the superconducting state belongs to one of the odd- 
parity representations of the tetragonal point group Z?4/j 12,15 . The quasi two-dimensionality of the Fermi liquid state 
in Sr 2 Ru04 makes it reasonable to assume that the momentum dependence of the superconducting state is described 
by functions f x (k) and f y (k) which obey the same symmetry transformation properties as k x and k y respectively 
under rotations of D^h (but otherwise are arbitrary). When there is no spin-orbit coupling, the spin-triplet state has 
six- fold degeneracy 20,21 , the transition temperature T c will be same for any linear combination for gap functions given 
in Table 1. When spin-orbit coupling is included, the degeneracy of T c will be lifted. The stable state will either have 
the d-vector along c-axis or be a linear combination of in-plane d- vectors. 

To quantify the role of spin-orbit coupling, we follow an approach developed by Sigrist et al. 22 . In particular, the 
effect of spin-orbit coupling is included through the magnetic susceptibility. The Hamiltonian with a general pairing 
interaction is 

H = Y. e k C ks C -ks + \Y. Yl V fc,fc'; S1S2S3S 4 C L 1 C -fc S2 C -fc' S3 C fc' S 4' (2) 
k,s k,k' S 1 S 2S 3 S4 

where e^, is electron band energy measured from the Fermi energy and cj^ , Cf ts are the fermion creation and annihilation 
operators. As a concrete model, we use an effective pairing interaction that is due to spin fluctuations 22 . However, 
the results that we require later depend solely upon the splitting of the six-fold degeneracy (this can be incorporated 
in a model independent way within the quasi-classical theory). The effective pairing interaction we use is 

where / is an coupling constant, and Xnv(k, k') is the static susceptibility. The phenomenological form of x^„(fc, k') 
in a material of tetragonal symmetry is 

/ 9i(f x f x + fyfy)+9 2 (f X fL-fyfy) 93(f x f y + fyfL) \ 

V 9z{f X f' x + fyf'y) ) 

where gi,g 2 ,g3 and g z are phenomenological parameters. The self consistency equation for d- vector with the above 
interaction can be solved to get the superconducting transition temperature, ksT c = l.lAuj c exp{— 16// 2 7V(0)Vr) for 
the different representations T (Vr corresponds to the interaction for the representation T). These are listed in Table 
I. In the limit g\ = g z and g 2 = 53 = there is no spin orbit coupling and all the representations will be degenerate. If 



Rep (r) 


Gap Function 


Interaction (Vr) 


A lu 


xfx + yf v 


g z - 2(g 2 + g 3 ) 


A 2u 


%fy ~ yfx 


g z +2(52 +33) 


Bin 


xfx - yf v 


gz - 2(g 2 - g 3 ) 


B 2u 


&fy + yfx 


gz + 2(g 2 - g 3 ) 


E u 


z(fx ± ify) 


2gi - gz 



TABLE I: Gap functions and interaction strengths for the different representations of D ih . 

there is cylindrical symmetry then g 2 = .93. Notice that this does not imply that that all the in-plane pairing states are 
degenerate. For this to occur, gi = 53 = 0. It is instructive to use results from recent microscopic calculations to gain 
an insight into the relative size of {gi, g 2 , 33, g z } 21 ' 23 - Both these papers reveal that deviation from the isotropic limit 
is small, since all the representations have very similar transition temperatures. The results of Ref. 23 correspond to 
the limit 33 = 0, gi > g Zl and \g\ — g z \ » \g 2 \- Based on these results we will assume that \g 2 \, {g^l << \g\ — g z \. The 
case gi > g z corresponds to the chiral pairing state 21,22 . While the case g z > gi and \g 2 \, \g 3 \ « \gi — g z \ corresponds 
to the nearly degenerate in-plane rf-vector. We will consider both these cases in the following. 
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IV. EILENBERGER EQUATIONS FOR THE 7 BAND 

An important aspect for understanding the superconducting state in Sr 2 Ru04 is the band structure. In particular, 
the states near the Fermi surface are derived from the Ru ti g orbitals. The degeneracy of these orbitals are split 
by the tetragonal crystal field into a xy orbital and the degenerate {xz,yz} orbitals 24-27 . These two sets of orbitals 
have a different parity under a mirror reflection through the basal plane. Consequently, to first approximation, the 
7 sheet of the Fermi surface is comprised of xy Wannier functions, while the {a, (3} sheets of the Fermi surface are 
comprised of {xz,yz} Wannier functions. This leads naturally to orbital dependent superconductivity 28 ' 29 ; a theory 
for the superconducting state has different gaps on the {a,f3} and 7 sheets of the Fermi surface. This theory has 
experimental support through specific heat measurements in magnetic fields 30 ' 31 . These measurements indicate that 
for strong fields applied in the basal plane, superconductivity in the {a, (3} bands is suppressed and the 7 band has 
the dominant superconducting gap. In addition to these measurements, recent theoretical calculations indicate that 
the ratio of the {a, (3} band gaps to that of the 7 band gap is 0.15 in the high field limit 19 . Therefore, we restrict 
the following microscopic theory in the high field regime to a single band theory for the 7 band. 

Now we will explain briefly the approximate analytic solution of the fundamental quasi-classical equations for a 
single band within weak-coupling superconductivity under magnetic fields. Our notation and formulation follows that 
of Rcf. 19. The solution can be obtained by the following approximations: (i) the spatial dependence of the internal 
magnetic field is averaged by B, (ii) the vortex lattice structure is expressed by the Abrikosov lattice and (iii) the 
diagonal elements of the Green's function are approximated by the spatial average. In general, a magnetic field will 
mix different representations of the Dah- Consequently, the (/-vector will be a linear combination of the functions 
listed in Table I, 

d(R,k) = di(fl)[i/ x (fc) +yf y (k)} + d 2 (R)[xf y (k) + yf x (k)} + d 3 (R)[xf x (k) - yf v (k)] 

+d i (R)[xf y (k) - yf x (k)\ + z[ Vx (R)f x (k) + Vy (R)f y (k)}, (5) 

the form of f x (k) and f y (k) is given in next subsection. The approximation (ii) amounts to taking each order 
parameter component in the lowest Landau level, so that d(R, k) = A<j) (R)d(k), where </> , the lowest Landau level, 
is given by 

MR) = E c ^ lVnV ' cx P[-((^ - A V)/A) 2 /2], (6) 

n 

where p n — 2im/ (3, (3 being the lattice constant in the y direction, A = (2\e\B)~ 1 ^ 2 is the magnetic length, and the 
coefficients c„, which determine the type of vortex lattice, are such that < \(f>o(R)\ 2 >= 1 , A is the magnitude of gap, 
and d(fc) defines the angular dependence of d- vector. We have taken anisotropy into account by writing x = x'/x 1 ^ 2 
and y = y'x 1 ^ 2 - For a conventional superconductor, even though all approximations mentioned above are valid near 
H c2 , comparisons with reliable numerical calculations suggest that the solution is competent quantitatively in wide 
region of the {T,H) phase diagram except in very low T and H regions 32 . For further detail we recommend the 
reader to refer to the literature 32-36 . In principle, we should consider higher Landau level solutions in this problem. 
However, within Ginzburg Landau theory our solution is exact in the high field limit provided the magnetic field is 
in the basal plane. This indicates that it is reasonable to keep only the lowest Landau level (for the field along the 
c-axis, other Landau levels must be included 13 ). The expression for the free energy measured relative to the normal 
state energy 32 , which is given for strongly type-II superconductors, B ~ iJ, in the clean limit is, 

fWiVo = In (^) (M 2 + M 2 ) + E I*! 2 In (|r) + 2^r£ - (/)) , (7) 

where N is the total density of states (DOS) in the normal state, and ui n = (2n + l)nT is the Fermionic Matsubara 
frequency. Eq. 7 generalizes the corresponding expressions in earlier works 32-36 by including more than one irreducible 
representation. The function / is given by 

1 + 9 \vj_(k) J 

with 

-1/2 

(9) 



(8) 



.9 = 



1 + 



2A 



A 2 |^(fc)| 2 W/V„), 



■5 



where u n — 2Au> n /v±(k), W(z) = e~ z crfc(— iz) is the Faddeeva function, and |^s(fc)| 2 = d(fe) • d*(fe). Here v±(k) 
is the component of v perpendicular to the field, which for an in-plane field of the form H = H(cos(j) h ,sin(f> h7 0), is 
given as 

5i(fe) = X - 1/2 v 2 z + x 1/2 v 2 F sin 2 (0 - fa), (fO) 

here \ = X v f/v c is an anisotropy parameter (we let \ be arbitrary). We have taken v z (k) — u c sgn(fc z ). All fields 
are measured in units 2n 2 T^ /(evFV c ) and h and ks are both 1. We will use this formalism to examine the in-plane 
magnetic field phase diagram for the chiral p-wave state and the nearly degenerate in-plane d- vector separately. 

A. Momentum dependence of the gap function 

To complete the description of the superconducting state, we must specify the functions fi(k) (i = {x,y}). The 
most general gap function consistent with translational invariance and with the appropriate rotational properties of 
D 4h is 

oo 

fi( k ) = c «,™ sin ( nfc ») cos(mfcj), i 7 j=x,y i^j, (ff) 

n— l,m— 

where c n m are complex coefficients, ki in units it /a, and a is the lattice spacing. Here n — f, m = represents 
a Cooper pair formed by nearest neighbor (nn) interactions and n = 1, m = 1 to a Cooper pair formed from next 
nearest neighbor (nnn) interactions. In general, increasing m, n corresponds to forming Cooper pairs from interactions 
between increasing number of neighbors. We will restrict ourselves here to nn and nnn pairing interactions. This has 
some support from microscopic calculations. In particular, the theory for nn interactions was originally proposed by 
Miyake and Narikiyo 37 and was also examined by Nomura and Yamada 38 . Cooper pairs for which nn interactions are 
not important but which have a substantial contribution from nnn interactions have been proposed recently by Arita 
et al. 39 with the assumption of large on-site and nearest neighbor Coulomb interactions. In particular, the latter paper 
proposes a gap function of the form d = z[sin(k x + k y ) — i sin(k x — k y )} = z\/2e _ "/ 4 [sin k x cos k y +i sin k y cos k x ]). We 
take Cooper pairs to be formed by nn and nnn interactions and take the form of f x to be, f x = sin(fca;)[l + ecos(k y )]. 
We keep the parameter e to be arbitrary and allow it to vary. We take the Fermi surface to be cylindrical with 
(k x , k y ) — ttR(cos (j>, sin <j>) where R = 0.9 approximates the Fermi surface of the gamma band (in the third subsection 
of next section we will take R=0.79 for reasons that will be apparent). 

V. CHIRAL p-WAVE STATE 

For the chiral p-wave state, we use d(fc, R) = zA(f> a (R)ip(k) with 

ip{k) = [cos9f x (k) + sm6e iC f y (k)], (12) 

(not normalized for notational simplicity). For this gap function, the Ginzburg Landau order parameter takes the 
specific form r)(R) = A</> (.R)(cos#, e 1 ^ sin 9). 

A. Upper critical field 

The anisotropy in upper critical field for (1,0,0) and (1,1,0) directions is a generic feature of the E u theory. Since 
experimentally it has been observed that the upper critical field is relatively isotropic for in-plane fields 14 , we ask if 
it is possible to reproduce this. The only free parameter in the theory is e which describes the anisotropy in the gap 
function. Surprisingly, we have found that it is possible for e = 10. The upper critical field for the field in the (1, 1, 0) 
and the (1, 0, 0) directions for the three values of e; e = 0, e = oo, and e = 10 for the chiral superconducting state has 
been shown in Figs. 1,2, and 3. In these figures we also show the stable solutions for the order parameter at H C 2. 
For a field along the (1, 0, 0) direction, the possible stable solutions are r\ = (1, 0) or r] — (0, 1). For a field along the 
(1, 1,0) direction, the possible stable solutions are r] — (1,1) or 77 = (1,-1). 
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FIG. 1: Upper critical fields for e = for the fields along the (1,0,0) and (1,1,0) directions. The field in units 2tv 2 T 2 /(evFV c ) 
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FIG. 2: Upper critical fields for e = oo and for the fields along the (1,0,0) and (1,1,0) directions. 



B. Phase diagram 



Another generic feature of the E u theory is that multiple vortex phases exist for in-plane fields. For a field along 
the (1,0,0) direction, the solution near H C 2 is r\ — (0,1) (for the three values of e discussed above this was the 
case) then as field is reduced for fixed temperature, a second transition occurs at H 2 for which the rj = (1,0) 
component becomes non-zero. Such phase transitions have only been examined within Ginzburg Landau theory 28,40 . 
The Eilenberger equations discussed above allow for the examination of this phase transition throughout the entire 
H-T phase diagram. Here, we apply this approach to the e = 10 gap function for a (1, 0, 0) field direction, the resulting 
phase diagram is shown in Fig. 4. 

The specific heat as a function of temperature is also shown in Fig. 5 for H/ H® 2 = 0.21. This plot clearly shows 
a second transition in the specific heat. Such a second transition has not been seen in specific heat measurements. 
This represents a difficulty for the E u theory. In general, we have not been able to find a microscopic theory that can 
account for both the lack of anisotropy in the upper critical field and the lack of the second transition. It is possible 
that experiments have not seen the predicted specific heat jumps due to the broadening associated with fluctuations 
in the vortex phase, or due to sample inhomogcneitics. In the case of UPt3, for which multiple phase transitions 
in the vortex phase have been observed, the entire phase diagram was found through ultrasound measurements 41 . 
Specific heat measurements mapped out some portions of the phase diagram 42 but they did not show clear anomalies 
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FIG. 3: Upper critical fields for e = 10 and for the fields along the (1,0,0) and (1,1,0) directions. The two upper critical fields 
almost identical. 
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FIG. 4: Phase boundaries for e = 10 and for the field along the (100) direction. 

throughout the entire phase diagram 3 . Therefore, it would be useful to look for such transitions in the phase diagram 
for Sr 2 RuC>4 with other probes such as ultrasound. 



C. Tetracritical point 

While there has been no evidence of multiple superconducting transitions in the low-field regime, two superconduct- 
ing transitions have been observed in the high- field range 14, . In one aspect these transitions are natural candidates 
for two transitions discussed above. In particular, the vanishing of the second transition (which does not occur at H c2 ) 
as the field is rotated away from the in-plane direction is consistent with the above predictions. However, the second 
transition is only observed for T/T c < 0.1 and appears to intersect the upper critical field line. This is inconsistent 
with the above prediction which predicts that this transition should exist for all temperatures. Here we explore a 
possible explanation for this transition that is based on results of the Eilenberger equations. 

We have found that for small parameter ranges in the model described above, it is possible that the solution for 
the order parameter at the upper critical field changes as a function of temperature. In particular for a field along 
the (1, 0, 0) direction, the low temperature solution (at H C 2) is rj = (1, 0) and there is a transition as temperature is 
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FIG. 5: Specific heat as a function of temperature for a fixed field H/H° 2 = 0.21 and for e = 10. 

reduced so that the solution at H c2 becomes rj = (0, 1). The phase diagram that emerges bears a striking similarity 
to the observed results. 

In Figs. 6 and 7, we show the phase diagram and specific heat calculations for R = 0.79 and e = —12. Note that a 
negative e corresponds to a repulsive interaction between Cooper pairs formed from nearest and Cooper pairs formed 
from next nearest neighbors. This theory would still require the existence of two transitions up to T = T c . However, 
as this phase diagram shows, the two transition lines between T c and the temperature of the tetracritical point are 
very close to each other and will be very difficult to observed experimentally. 




T/T c 



FIG. 6: Phase diagram showing a tetracritical point. Between T c and the tetracritical temperature there are two phase 
transitions, shown in the inset, that are difficult to distinguish from each other. 

While this phase diagram agrees with that observed experimentally for a field along (1,0,0) direction, this choice of 
parameters also exhibits a moderate w 15% anisotropy of the in-plane upper critical field. Furthermore, for the field 
along the (1, 1,0) direction the phase diagram resembles that of Fig. 4 (there is no tetracritical point). Therefore, 
these results can only be taken as suggestive since this set of parameters cannot account for all observed features. It 
is possible that the gap on the {a, f3} bands may improve the agreement between theory and experiment. It appears 
that at low fields these gaps cannot be neglected 19 and the suppression of these gaps relative to that of the 7 may 
provide a more robust mechanism for the appearance of a tetracritical point in the relatively low field regime. This 
can occur if these bands prefer orthogonal order parameter solutions at H c2 , as was often the case for calculations for 
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FIG. 7: Specific heat as a function of temperature for fixed magnetic fields. 



different E u gap functions in Rcf. 17. A second possible explanation for both the observed lack of anisotropy and the 
existence of a tetracritical point is that the parameter e may not be constant for all temperatures and magnetic fields 
as we have used here. A complete description would require an effective two-gap theory for which e is determined 
self-consistently. In such a theory e can develop a temperature and magnetic field dependence (though in many 
circumstances, e will be approximately constant). A complete analysis of this would require a detailed knowledge of 
the microscopic interactions giving rise to superconductivity. This is beyond the scope of this paper. 



The chiral p-wave state has difficulties explaining both the observed absence of anisotropy in H c2 and the absence 
of additional phase transitions in the vortex phase. For this reason we also consider the nearly degenerate in-plane 
d-vector. Initially we consider under what circumstances the nearly degenerate d-vector is consistent with other 
experiments (Knight shift, /iSR, vortex lattice structure, and Josephson experiments). Then we examine the specific 
heat as a function of magnetic field in the vortex state and show that the second anomaly is rapidly suppressed by 
the breaking of the degeneracy of the four components of the in-plane d-vector. 

Knight shift measurements can be naturally accounted for by an in-plane d-vector. The most recent observation 
is that for magnetic fields along c-axis there is no change in the spin susceptibility as temperature is reduced 6 . The 
simplest interpretation of this result is that the d-vector is in the basal plane. The earlier observation that the spin 
susceptibility is unchanged for in-plane magnetic fields would require that the d-vector is free to rotate in the basal 
plane. This implies that all the in-plane d-vector states are degenerate or nearly degenerate. The observed square 
vortex lattice for the field along the c-axis is also consistent with a degenerate in-plane d-vector. This will follows 
from an Ginzburg Landau analysis, where it can be shown that free energy for the degenerate in-plane d-vector has 
equilibrium solutions with the same properties as those of the Ginzburg Landau theory for the chiral p-wave state. 

The muon spin relaxation (/iSR) measurements of Luke et al. 4 have found an increased spin relaxation rate in 
the superconducting state with zero applied magnetic field. This has commonly been interpreted as evidence for a 
superconductor that breaks time reversal symmetry in the bulk. However, any bulk internal magnetic field must be 
screened due to the Meissner effect. Consequently, /xSR only probes internal magnetic fields due to inhomogeneities 
such as impurities or domain walls between degenerate superconducting states. It has been shown that a supercon- 
ducting state that does not break time reversal symmetry in the bulk can still give rise to local internal fields 15 . The 
important condition for such internal magnetic fields to exist is that the superconducting order parameter has more 
than one degree of freedom. This distinction is emphasized here because for nearly degenerate in-plane d-vectors, the 
bulk superconducting state in zero applied field does not break time reversal symmetry. This does not imply that 
such a state is inconsistent with /xSR measurements. However, it does require that the different in-plane d-vector 
representations are nearly degenerate (the different T c values must lie close to each other). 

The most difficult experiments to explain with an in-plane d-vector are the Josephson experiments. The most 
recent of these has found that for Sr 2 Ru04-Auo.5ln .5 SQUID, there is a ir phase difference in the Josephson current 
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when the two junctions have opposite normals 8 . While this is generally expected for a p-wave superconductor, it 
cannot be explained by an in-plane d-vector. Such a d-vector does not allow for a Josephson current between an 
odd-parity superconductor and an isotropic (s-wave) superconductor when the junction has a normal perpendicular 
to the Sr 2 Ru04 c-axis. The existence of such a Josephson current implies a d-vector aligned along the c-axis. An 
explanation of such a Josephson current within an in-plane d-vector approach would therefore require that at an 
interface, the d-vector is along the c-axis and in the bulk it is in-plane. A similar scenario has been proposed by 
Bahcall in the context of the cuprate superconductors (in this case, the order parameter near the interface is s-wave 
and becomes cZ-wave in the bulk) 45 . In support of such a picture, the d-vector is almost certainly along the c-axis for 
an interface with a normal perpendicular to the c-axis. This is a natural consequence of a stronger spin-orbit coupling 
at the interface than in the bulk (spin-orbit coupling is governed by the gradient of the single-particle potential). 
It is the spin-orbit coupling that governs the orientation of the d-vector. A Rashba spin-orbit coupling of the form 
ajifi ■ k x S(k) (where h is the interface normal, k is the fermion wave-number, S(k) is the fcrmion spin, and an is 
a coupling constant) would give a d with a component along the c-axis if n lies perpendicular to the c-axis 46 . If the 
bulk d-vector is in-plane and the d-vector lies along the c-axis near the interface, then an analysis following that of 
Bahcall would imply that the 7r squid experiment of Nelson et al. should sometimes see a ir phase shift and sometimes 
no phase shift. Nelson's data indicates that there is always a n phase shift. However, given that data on only three 
samples are presented, it may be prudent to await further results before ruling out an in-plane d-vector in the bulk 
on the basis of these experiments. 

The general form for an in-plane d-vector will be a linear combination of four different in-plane representations 
listed in Table I, 

d(R,k) = d!{R)[xf x {k) + yf v (k)} + d 2 (R)[xf y (k) + yf x (k)} + d 3 (R)[xf x (k) - yf y (k)} 

+d 4 (R)[xf y (k)-yf x (k)}. (13) 

We parameterize [di(R), d 2 (R), ds(R), d4(R)] as A(f>o(R)(cosip cos9,simpsm(f> 7 cos ^ sin ^, sin ^ cos 0) so that d- 
vector can be written as 

d(R,k) = A(R)[x<p x (k)+y<p y (k)}, (14) 

with (p x (k) = cos i(j(cos 9 + sin9)f x (k) + sin i(j (sin <p + cos <j))f y (k) and <p y (k) = cosf/^cosfl — sin 9) f y (k) + sin?/;(sin</> — 
cos<p)f x (k). With these definitions \f(k)\ 2 = \(p x (k)\ 2 + \(p y (k)\ 2 . 

Prior to presenting the results for the nearly-degenerate in-plane d-vector, we briefly give the results for the 
degenerate in-plane d-vector. This corresponds to the situation that g 2 = 53 = in Table I. In zero applied field 
there are many degenerate solutions for this case. In particular, all the solutions in Table I and the solutions 
d(k) = e(k x ± ik y ) 7 where e is any unit vector in the basal plane, are all degenerate ground state solutions 20 . If 
a magnetic field is applied along any of the two- fold symmetry axes, then there will be two transitions as field is 
reduced (as there was in the chiral p-wave case) . Unlike the chiral p-wave case the second transition can occur in two 
ways. To illustrate this, consider an applied field along (0, 1,0) direction, the high field state (the state corresponding 
to the transition from normal state to superconducting state) will be either xk x or xk y (it has been assumed here 
that the d-vector prefers to be perpendicular to the magnetic field). Consider xk y to be concrete. The second 
transition will appear as magnetic field is reduced for fixed temperature. The second transition exists because of 
the appearance of either a xk x component [the corresponding zero-field ground state will be x(k x ± ik y )]; or a yk x 
component (the corresponding zero-field ground state will then be xk y ± yk x ). Strictly speaking, the latter transition 
will be energetically less favorable because the d-vector is not perpendicular to the magnetic field. However, it is 
the latter transition that will play a more important role when the degeneracy between the four in-plane d- vectors is 
broken. In this case the solutions at zero-field belong to a single irreducible representation while the other zero-field 
solutions x(k x ± ik y ) belong to a mixture of more than one irreducible representation. 

In Fig. 8, we show the specific heat as a function of temperature for different values of g 2 (we have set g 2 = 33 in 
the following) with fixed magnetic field, H/H® 2 = 0.21 along (0,1,0) direction for e = 10 and R = 0.9. For g 2 = 0, 
as discussed above, the second transition will exist and there are two specific heat anomalies. The second transition 
is removed by a finite value of g 2 /g z - The key result is that the anomaly for the second transition is very quickly 
suppressed by a non-zero g 2 /g z - Note that the anisotropy in H c2 will still be small for small values of g 2 . Consequently, 
a nearly degenerate in-plane d-vector can explain the existing experimental observations on the H-T phase diagram. 

The nearly degenerate in-plane d-vector can account for the in-plane phase diagram and can qualitatively account 
for other key experimental results in Sr2Ru04. However, prior to carrying out further calculations with this state 
we note that it should be possible to rule such a state out experimentally in the near future. In particular, there 
are two predictions that can be made about an in-plane d-vector. The first has been mentioned above: further n 
SQUID experiments should reveal the existence of squids with no phase shift as well as squids with ir phase shifts. 
Also, further Knight shift experiments should show a suppression in the spin susceptibility for low enough in-plane 
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FIG. 8: Specific heat for e — 10 for the field along the (1,0,0) direction, with fixed H/H^ 2 = 0.21 for different values of g2 
(which is measured here in units g z ). 

magnetic fields. This will occur because in zero-field the d-vector will correspond to a single-component irreducible 
representation once gi and t/3 are non-zero and therefore contain a component that is along the applied field. 

VII. CONCLUSIONS 

To address an apparent conflict between theoretical predictions of chiral p-wave (the E u representation) theories of 
the superconducting state in Sr 2 Ru04 and the lack of corresponding observations, we have carried out quasi-classical 
calculations of the superconducting phase diagram for in-plane magnetic fields. This has been done for both the 
chiral p-wave state and for the nearly degenerate in-plane d-vector. For a gap function with momentum dependence 
due to a combination of nearest and next nearest neighbor interactions defined on the 7 band, we find that a small 
anisotropy in the upper critical field as the field is rotated in plane is possible. However, the same gap functions give 
rise to an additional phase transition in the vortex state which has not been observed experimentally. For a narrow 
range of parameters, the theory gives rise to a tetracritical point in the H-T phase diagram. When this tetracritical 
point exists, the resulting phase diagram closely resembles the experimentally measured phase diagram for which two 
transitions are only observed in the high field regime. We have also argued that an in-plane d-vector that can easily 
rotate in the basal plane is consistent with existing experimental results. 
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